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STATISTICAL  ANALYSIS  OF  KEVLAR  A9/EPOXY  COMPOSITE  STRESS-RUPTURE  DATA* 


Abstract 

Statistical  analyses  are  presented  for  LLNL  stress-rupture  data  sets 
involving  Kevlar  49/epoxy  strands  and  NASA  Kevlar  49/epoxy  spherical  pressure 
vessels  subjected  to  sustained  loading.  Raw  data,  summarized  inferences,  and 
figures  are  included. 

Introduction 

Lifetime  properties  of  Kevlar  49/epoxy  composites  subjected  to  sustained 
loading  may  be  extracted  from  the  statistical  analysis  of  stress-rupture  data 
on  composite  strands  and  NASA  spherical  pressure  vessels  generated  at  LLNL  in 
the  1970s.  The  purpose  of  this  report  is  to  present  the  collection  of 
statistical  analyses  made  on  these  data  using  the  Weibull  maximum  likelihood 
accelerated  testing  methology  described  in  Reference  4. 

The  following  stress-rupture  data  sets  are  analyzed  here: 

0  Kevlar  49/epoxy  strands  at  room  temperature  (UV  light  present) 

0  Kevlar  49/epoxy  strands  at  elevated  temperatures  (no  UV  present) 

0  Kevlar  49/epoxy  NASA  spherical  vessels  at  room  temperature  (no  UV 
present) 

Included  in  the  report  is  a  complete  collection  of  available  relevant  raw 
data,  i.e.,  static  strength  measurements  and  lifetime  measurements  (exact, 
censored,  and  grouped),  as  well  as  summary  descriptions  of  the  materials, 
e.g.,  mean  fiber  weight,  cross-sectional  area,  fiber  volume  content,  etc. 
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A.  Kevlar/epoxy  strand  data  analysis:  lifetime  versus  stress  (room 


temperature  environment  with  UV  light). 

1 .  Materials 

The  Kevlar  49/epoxy  strands  used  in  room  temperature  strength  and 
lifetime  tests  are  described  in  Table  1.  Added  detail  and  descriptions  of  the 
testing  equipment  may  be  found  in  References  5  and  7. 


Table  1.  Summary  of  material  characteristics. 
(Room  temperature  strands) 


construction:  filament-wound  from  single  end  (380  denier,  approximately  267 
filaments)  pre-production  Kevlar  49  fiber  without  finish 


epoxy  matrix 

0  Union  Carbide  ERL  2258/ZZL  0820 
0  weight  ratio  100/29 

0  vacuum  impregnation  of  fiber  with  epoxy 


fiber  specific  gravity: 
mean  fiber  weight: 
mean  cross-sectional  area: 
curing: 


1,45 

0.0409  g/m  (based  on  5  specimens) 

-5  2 

4.364  X  10  in  (based  on  5  specimens) 
3  h  at  93‘’C,  2  h  at  163®C 


mean  fiber  volume  content:  71.5%  (based  on  35  specimens) 

length  of  specimens  for  lifetime  and  strength  tests:  10  in 
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2.  Data 


A  strength  test  was  performed  on  53  specimens  resulting  in  an 
estimated  ultimate  tensile  strength  of  505.7  ksi,  with  estimated  coefficient 
of  variation  2.468%.  The  complete  set  of  breaking  loads  is  presented  in 
Appendix  A.  Stress-rupture  (lifetime)  tests  were  subsequently  carried  out  at 
room  temperature  in  the  presence  of  UV  light  over  an  eight  year  period  on  a 
total  of  553  strands  at  seven  stress  levels  ranging  from  90%  down  to  50%  UTS. 
The  data  are  summarized  in  Table  2. 

Table  2.  Room  temperature  strand  data  summary. 


stress,  ksi 

(%UTS) 

#strands 

#exact 

failure 

times 

^grouped 

failure 

times 

#censored 

failure 

times 

marginal 

A 

a 

ML  est. 

A 

b 

455.2 

(90.0) 

101 

99 

2 

0 

1.09 

-.0159 

440.0 

(87.0) 

100 

100 

0 

0 

1.04 

1.37 

424.8 

(84.0) 

103 

103 

0 

0 

1.09 

3.09 

404.6 

(80.0) 

100 

100 

0 

0 

.923 

5.37 

354.0 

(70.0) 

49 

47 

2 

0 

.496 

9.20 

303.5 

(60.0) 

50 

41 

9 

0 

.292 

10.7 

252.9 

(50.0) 

50 

4 

0 

46 

.458 

12.3 

A  marginal  ML  estimate  for  a  given  stress  level  is  an  estimate  whose 
computation  is  based  solely  on  the  data  from  that  stress  level.  An  exact 
failure  time  is  a  lifetime  recorded  as  a  point,  e.g.  1456.3  hours.  A  grouped 
failure  time  is  a  lifetime  recorded  as  an  interval,  e.g.  [1450,  1474]  hours. 
Such  is  necessitated  when  the  precise  failure  time  of  a  specimen  is  not 
recorded  exactly,  although  the  time  is  known,  as  in  a  timer  breakdown,  to  be  a 
number  between  two  specified  limits.  If  a  non-failed  specimen  is 
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removed  from  test  independently  of  its  condition,  say  because  the  strand  slips 
from  its  clamp,  or  an  earthquake  destroys  the  experiment,  or  the  decision  is 
made  to  terminate  all  testing,  the  removal  time,  called  censoring  time,  is 
recorded . 

In  the  next  section  the  collection  of  exact,  grouped,  and  censored 
failure  times  with  accompanying  stresses  is  analyzed  statistically  to 
characterize  the  lifetime  distribution  of  strands  as  a  function  of  stress. 

3.  Statistical  Analysis  of  the  Data 

The  seven  Weibull  probability  paper  plots  of  the  lifetime  data 
partitioned  according  to  the  individual  stress  levels  support  the  basic 
assumptions  made  here 

(i)  that  the  lifetime  distribution  for  a  given  stress  level  a  is 
two-parameter  Weibull,  and 

(ii)  that  the  Weibull  parameters,  a  (shape)  and  8  (scale, 
characteristic  life),  depend  on  a. 

The  statistical  methodology  introduced  and  used  in  the  analysis  of  S-glass/ 
epoxy  strand  stress-rupture  data  (Reference  2)  is  accordingly  appropriate  also 
in  the  Kevlar  setting.  Thus,  the  lifetime  T^  of  a  Kevlar  49/epoxy  strand 
under  constant  stress  a  is  modeled  by  the  probability  density  function, 

fj  (t)  =  [a(a)/6(a)][t/e(a)]“^^^”^exp{-[t/6(a)]°‘^^^/,  t  >0  . 

a 

Subsequent  computations  are  simplified  by  consideration  of  Y  =  fin  T  ,  whose 

a  a’ 

probability  distribution  is  that  of  a(a)W  +  b(a),  where  a(a)  =  l/a(o), 
b(a)  =  fin  8(o),  and  W  has  the  extreme  value  density. 
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=  exp  (w  -  e'^)  ,  -  «>  <  w  <  “  . 

The  effect  of  the  transformation  =  £n  is  a  reparametrization  of 

a(o)  and  3((J)  to  a(a)  and  b(a)  which  behave  as  extreme  value  scale  and  location 

parameters,  respectively. 

Examination  of  marginal  maximum  likelihood  (ML)  estimates  of  the 
parameters  (a,  b)  for  the  seven  individual  stress  levels  (see  Table  2) 
suggests  that  the  Weibull  parameter  to  stress  dependency  (ii)  can  be  modeled 
adequately  by  the  polynomials 

a  (a)  =  6^ 

b(a)  =  6^  +  65O  +  6^0^  +  . 

Properties  of  the  lifetime  distribution  for  a  given  stress  a  are  thereby 
expressible  in  terms  of  the  polynomial  coefficients  6^^,  0y.  Moreover, 

estimates  of  lifetime  distribution  characteristics  of  interest  may  be  generated 
from  estimates  of  6^,  ...,  6^. 

Maximum  likelihood  estimates  of  the  coefficients  9^^,  ...»  6y  were 
computed  from  the  aggregated  collection  of  553  exact,  grouped,  and  censored 
failure  times  recorded  in  terms  of  hours  for  the  seven  stress  levels  measured 
in  ksi.  The  results  obtained  are 

(1)  (6^,  02,  Qy  e^,  6y  e^,  6^)  =  (1.22750,  -  9.68792  x  lO"^,  2.16947  x  10 

3.23296  X  10^,  -2.18769  x  lO"^, 

8.11080  X  10”'^,  -1.07399  x  lO"^) . 
Accordingly,  the  ML  estimates  of  the  functions  a (a)  and  b(a)  are 
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(2)  a(a)  =  0^  +  e^cT  +  and  b(a)  =  0^  +  , 

and,  equivalently,  the  ML  estimates  of  the  dependencies  a{o)  and  B(a)  are 
a(a)  =  l/a(a)  and  6(a)  =  exp[b(a)].  Of  particular  interest  are  the 
quantities  tp(a)  and  R|.(a),  respectively  the  p^*^  quantile  and  the 
reliability  at  time  t  of  the  lifetime  distribution  for  stress  a;  t  (a)  is 
the  time  corresponding  to  a  failure  probability  of  p  (a  strand  at  stress  a  will 
fail  before  time  tp(a)  with  probability  p),  and  R^(a)  is  the  probability  of 
surviving  past  time  t  (the  lifetime  of  a  strand  at  stress  a  will  exceed  t  with 
probability  R^(a).  It  is  seen  that  tp(a)  =  exp[yp(a)],  where  yp(a)  = 
s(a)  Wp  +  b(a),  and  Wp  =  jin(-jln(l-p)).  Also  R^(a)  =  exp{-exp[(s,n  t  - 
b(a))/a(a)]}.  The  ML  estimates  of  tp(a)  and  R^(a)  are  therefore 

(3)  tp(a)  =  exp[yp(a)]  ,  where  yp(a)  =  a(a)  Wp  +  b(a),  and 

(4)  R^(a)  =  exp{-exp[(j,n  t  -  b(a))/a(a)]}  . 

The  estimated  standard  errors  (estimated  standard  deviations)  of  all  ML 
estimates  given  here  are  computable  from  results  presented  in  Appendix  A. 

The  maximum  likelihood  estimates  of  selected  quantiles  and 
reliabilities  as  functions  of  stress  are  displayed  in  the  accompanying 
Figures  1  through  5.  In  Figure  1,  the  data  are  presented  in  raw  fashion  as 
sample  quantiles.  At  each  of  the  seven  experimental  stress  levels,  base  10 
logarithms  of  the  standard  nonparametric  estimates  of  selected  quantiles  are 
plotted.  For  a  given  quantile  probability  p  and  experimental  stress  level 
Oi  this  so-called  Kaplan-Meier  estimate  is  essentially  the  time 
corresponding  to  a  failure  proportion  of  p  in  the  sample  of  failure  times  at 
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level  a.  For  the  selected  quantile  probabilities  p  =  .02,  .04,  .05,  .10, 

.30,  .50,  and  .90,  these  estimates,  where  computable,  are  connected  by  line 
segments.  The  result  is  a  crude  depiction  of  the  relationship  between  stress 
and  failure  time.  In  Figure  2,  base  10  logarithms  of  the  maximum  likelihood 
estimates  of  the  same  selected  quantiles  are  plotted  against  stress  (i.e.,  log 

A 

tp(a)  vs.  a).  The  general  agreement  of  the  nonparametric  and  ML  plots 

is  shown  in  Figure  3,  where  the  competing  estimates  are  superimposed  for  the 

quantile  probabilities  .02,  .05,,  .10,  and  .50. 

Figures  4  and  5  give  plots  of  estimated  quantiles  and  reliabilities 

for  situations  of  potential  engineering  use.  In  Figure  4,  base  10  logarithms 

of  the  ML  estimates  tp(a)  are  plotted  for  the  quantile  probabilities 

p  =  10  10  10  10"^,  10"^,  10"L  and  .50  in  the  range  100  to  450  ksi 

(19.8%  to  89.0%  UTS).  In  Figure  5,  base  10  logarithms  of  the  ML  estimated 

failure  probabilities,  1  -  R^(a),  are  plotted  for  the  times  1,  3,  5,  10,  15, 

20,  25,  and  30  years  in  the  range  125  to  275  ksi  (24.7%  to  54.4%  UTS).  To 

illustrate  the  use  of  these  figures,  consider  estimating  the  p  =  .01-quantile 
s  t 

(1  percentile)  of  the  distribution  of  lifetimes  of  strands  subjected  to 
constant  stress  a  =  200  ksi.*  From  Figure  4,  log  t  qi(200)  =5.1  so  that 
t. 01 (200)  =  10^*  =  126,000  hours  =  14. 4  years.  For  greater  accuracy,  the 

ML  estimate  can  be  computed  exactly  from  the  equations  (1),  (2),  and  (3).  Here 
y  q^(200)  =  11.7016  so  that  t^Q^(200)  =  exp(11.7016)  =  120,765  hours  =  13.8 
years.  From  results  presented  in  Appendix  A,  uncertainty 


The  inference  applies  only  to  the  population  of  strands  equivalent  in 
fabrication  to  those  tested  and  subject  to  an  environment  equivalent  to 
that  controlled  in  the  experiment. 
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in  the  statistical  estimation  procedure  can  be  quantified  to  give  an 
approximate  90%  confidence  interval  for  t  gj^(200)  of  [6.82,  27.8]  years. 

Now  consider  estimating  the  reliability  at  10  years  for  strands  at  stress  a 
=  200  ksi.*  From  Figure  5,  log  (1  -  y^(200))  i  -2.95,  which  implies 

^10  yr^^*^^^  ^  .9989.  The  exact  value,  computed  from  equations  (1),  (2),  and 
(A),  is  .9987.  The  corresponding  approximate  90%  lower  confidence  bound  as 
described  in  Appendix  A  is  computed  to  be  .9926. 

A.  Discussion 

The  lifetime  distribution  of  Kevlar  A9/epoxy  strands  as  a  function 
of  stress  has  been' estimated  by  maximum  likelihood  methods.  Formulas  (1)  - 
(A)  and  Figures  1,  A,  and  5  provide  ML  estimates  of  various  quantiles  and 
reliabilities  for  a  large  range  of  stress  levels.  Standard  errors  of  all  ML 
estimates  are  computable  from  results  given  in  Appendix  A.  The  usefulness  and 
applicability  of  these  statistical  inferences,  however,  are  limited  by  the 
following  sources  of  uncertainty. 

(a)  Sampling.  Because  lifetimes  recorded  at  the  given  experimental 
stress  levels  obey  probability  distributions,  the  actual  data  aggregated  will 
be  diverse  and,  especially  for  small  sample  sizes,  not  necessarily 
representative  of  the  underlying  populations.  The  degree  of  uncertainty 
attributable  to  sampling  variation  for  the  ML  estimates  generated,  although 
quantified  by  the  appropriate  estimated  standard  errors,  can  be 
unsatisfactorily  great  (i)  in  cases  of  extrapolation  in  stress  level  far 
beyond  the  data,  and  (ii)  in  estimation  of  extreme  tail  distribution 
parameters,  e.g.,  or  t2g-6. 


(b)  Weibull  distribution  assumption.  The  Weibull  model  is  an  ideal 
that  is  only  approximately  realized  in  experimentation.  To  the  degree  the 
Weibull  assumption  is  violated,  ML  estimates  are  biased  away  from  the  unknown 
parameters  being  estimated.  Such  bias  is  more  substantial  for  extreme  tail 
parameters,  e.g.  t^Q-3,  or  t-jg-e,  than  for  the  median,  t  g.  Standard  errors 
of  all  ML  estimates  also  are  erroneous  to  the  extent  the  Weibull  assumption  is 
violated. 

(c)  Assumed  relationship  of  Weibull  parameters  to  stress.  The 
assumed  stress  models  which  specify  the  functional  relationships  a(a)  and 
b(a)  have  a  crucial  bearing  on  all  ML  estimates.  The  selection  of  quadratic 
polynomial  and  cubic  polynomial  functions,  respectively,  for  the  LLNL  data 
provided  good  agreement  of  ML  estimates  with  standard  nonparametric  estimates 
(Figure  3)  and  the  highest  maximized  likelihood  among  all  functional 
relationships  considered.  In  short,  this  selection  fit  the  observed  data 
quite  well.  Unfortunately,  reasonableness  within  the  region  of  the 
experimental  stress  levels  does  not  guarantee  reasonableness  at  extrapolated 
stress  levels  far  beyond  this  region.  For  example,  two  sets  of  assumed 
functional  dependencies,  say  (a(a),  b(a))  and  (a*(a),  b*(a)),  may  give  rise  to 
respective  ML  estimates  of  parameters  which  within  the  experimental  stress 
region  are  nearly  coincidental.  Thus  t  Qi{a)  =  t  Qi(a)  for  all  a  between  250 
and  450  ksi.  However,  it  may  well  be  the  case  that  t  qi(100)  is  vastly 

A  ■jUi’ 

different  from  t  g^(lOO),  and  the  true  value,  t  gi(lOO),  may  be  far  away  from 
either  estimate.  Misspecif ication  of  the  functions  a{cr)  and  b(CT),  like 
violation  of  the  Weibull  distribution  assumption  (b),  introduces  bias  into  the 
ML  estimates  and  their  standard  errors.  This  bias  may  be  arbitrarily  great  at 
extreme  stress  levels. 
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(d)  Extension  to  other  populations.  The  statistical  inferences 
based  upon  the  given  data  set  are  strictly  relevant  to  a  hypothetical 
population  of  strands  whose  characterization  in  terms  of  fabrication  and 
controlled  environment  is  equivalent  to  that  of  the  sampled  strands. 

Application  of  the  results  to  strands  made  from  a  different  process  or  epoxy 
formulation,  or  subjected  to  environments  with  different  temperature ,  humidity , 
or  light  conditions,  introduces  elements  of  uncertainty  which  at  this  point  are 
not  understood  well  enough  to  be  quantified. 

B.  Kevlar  49/epoxy  strand  data  analysis:  lifetime  versus  stress  and 
temperature  (environment  without  UV  light). 

1.  Materials 

The  Kevlar  49/epoxy  strands  used  in  elevated  temperature  lifetime 
tests  are  described  in  Table  3.  Added  detail  and  descriptions  of  the  testing 
equipment  may  be  found  in  References  6  and  7. 

2.  Data 

Kevlar  49/epoxy  strands  of  reportedly  the  same  batch  as  in  part  A 
were  subjected  to  stress-rupture  tests  at  various  elevated  temperatures  and 
selected  stresses.  No  light  was  present  during  testing.  The  combined  effect 
of  stress  and  temperature  on  lifetime  can  be  estimated  from  a  statistical 
analysis  of  the  observed  failure  times.  Since  the  experimental  elevated 
temperatures  fell  within  the  range  100-120°C,  it  was  decided  for  interpolative 
purposes  to  include  in  the  analysis  those  portions  of  the  room  temperature 
data  of  part  A 
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Table  3.  Summary  of  material  characteristics.  (Elevated  temperature  strands) 


construction:  filament-wound  from  single  end  (380  denier,  approximately  267 
filaments)  pre-production  Kevlar  49  fiber  without  finish 


epoxy  matrix 

0  Union  Carbide  ERL  2258/ZZL  0820 
0  weight  ratio  100/29 

0  vacuum  impregation  of  fiber  with  epoxy 


fiber  specific  gravity: 
mean  fiber  weight: 
mean  cross-sectional  area: 
curing: 

mean  fiber  volume  content: 
length  of  specimens  for  li 


1.45 

0.0432  gm/m  (based  on  10  specimens) 

-5  2 

4.619  X  10  in  (based  on  10  specimens) 

3  h  at  100“C,  2  h  at  170°C 
67.8%  (based  on  40  specimens) 
fetime  and  strength  tests:  225  irni  (8.9  in) 


which  were  felt  to  be  generally  unaffected  by  the  presence  of  UV  light,  namely 
the  four  highest  stresses.  (At  the  lowest  of  these  stresses  (80%  UTS),  all 
failure  times  were  less  than  1000  hours;  at  the  highest  stress  excluded  (70% 
UTS),  failure  times  were  as  high  as  two  years.)  Although  the  composite 
materials  under  both  elevated  and  room  temperature  conditions  had  essentially 
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identical  descriptions  (in  terms  of  epoxy,  denier,  curing,  fiber  volume, 
etc.),  measured  cross  sectional  areas  and  ultimate  tensile  strengths  differed 
slightly.  Accordingly  for  the  sake  of  compatibility  of  the  two  data  sets, 
nominal  experimental  stress  levels  were  given  minor  adjustments  to  achieve 
common  proportionality  to  actual  applied  loads,  and  an  overall  UTS  figure  of 
500  ksi  was  adopted.  The  combined  data,  consisting  of  1190  specimen  lifetimes 
at  a  total  of  17  stress/temperature  settings  are  summarized  in  Table  4. 

3.  Statistical  Analysis  of  the  Data 

Weibull  probability  paper  plots  of  failure  time  data  at  individual 
stress/temperature  settings  suggest  that  strand  lifetimes  may  be  modeled  well 
as  two-parameter  Weibull  random  variables  whose  parameters  a  and  g  depend 
on  both  stress  o  and  temperature  x.  As  in  the  analysis  of  part  A,  use  of  a 
logarithmic  transformation  achieves  the  parametrization  a(a,T)  =  l/a(a,T) 
and  b(a,T)  =  Jin  g{a,T).  The  stress /temperature  influence  on  a  and  b,  by 
examination  of  maximum  likelihood  estimates  of  a  and  b  for  the  17  individual 
experimental  settings  (see  Table  4),  is  seen  to  be  modeled  adequately  as 

a(a,T)=6^+0^+9  3T^  ,  and 

2  -1  -2  -1 

b(a,T)=04  +  9g7+0g:f  +07T  ■*’9  3^  ’ 

where  o  is  measured  in  ksi  units,  and  x  in  degrees  Kelvin  (273  +  degrees 
Centigrade).  The  interaction  term,  0^x'\  allows  the  incremental  stress 
effect  on  lifetime  to  be  different  at  different  temperature  levels  (and  also 
allows  the  incremental  temperature  effect  on  lifetime  to  be  different  at 
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Table  4.  Combined  elevated  temperature/room  temperature  data  summary. 


stress, 

ksi 

(%UTS) 

temperature 

°C 

#strands 

436.7 

(87.3) 

120 

48 

436.7 

(87.3) 

no 

39 

436.7 

(87.3) 

100 

46 

410.7 

(82.1) 

120 

46 

410.7 

(82.1) 

no 

78 

410.7 

(82.1) 

100 

53 

373.9 

(74.8) 

120 

80 

373.9 

(74.8) 

no 

76 

373.9 

(74.8) 

100 

76 

343.2 

(68.6) 

120 

73 

343.2 

(68.6) 

no 

54 

343.2 

(68.6) 

100 

58 

288.1 

(57.6) 

no 

59 

447.6 

(89.5) 

25 

101 

432.7 

(86.5) 

25 

100 

417.8 

(83.6) 

25 

103 

397.9 

(79.6) 

25 

100 

lexact  Igrouped  marginal  ML  est. 
failure  failure 


times 

times 

a 

b 

48 

0 

.670 

-3.27 

39 

0 

.476 

-2.26 

46 

0 

.482 

-1.77 

46 

0 

.568 

-1.71 

78 

0 

.477 

-1.25 

53 

0 

.738 

1 

• 

PO 

00 

o 

80 

0 

.607 

.364 

76 

0 

.754 

.757 

76 

0 

.915 

1.57 

73 

0 

.533 

2.08 

54 

0 

.592 

3.14 

58 

0 

.657 

4.00 

59 

0 

.446 

5.65 

99 

2 

1.09 

-.0159 

100 

0 

1.04 

1.37 

103 

0 

1.09 

3.09 

100 

0 

.923 

5.37 

14 


different  stress  levels).  The  model  for  the  transformed  shape  parameter,  a, 
is  simpler  here  than  in  the  previous  analysis  A,  a  possible  consequence  of  the 
absence  of  UV  light. 

The  maximum  likelihood  estimates  of  the  coefficients  e-j,  ...»  ©g, 
computed  (in  a  manner  analogous  to  that  of  analysis  A)  from  the  collection  of 
1090  failure  times  aggregated  from  the  17  stress/temperature  levels,  are  the 
fol lowing: 

(0|,  02,  03.  04,  05,  05,  07,  03,  0g)  =  (-7.94341  X  10"  ,  2.17931  x  10“^ 

5.23962  X  10^,  -1.02439  x  10^,  1.29433  x  10"\  -3.47076  x  10"^  , 

5.70114  X  10^  -4.21914  x  10®,  -6.15226  x  10^)  . 


The  ML  estimates  of  the  functions  a,  b,  a,  and  6  are  therefore  a(a,T)  =  ©i 
+  ©20+0  3T  »b(a,T)-04  +  0^+©^  +07T  ■*■08'^  > 

A  A  A  A 

a(a,T)  =  l/a(a,T),  and  3(a,T)  =  exp[b(a,T)].  Similarly,  ML  estimates  of  the 
lifetime  parameters  tp(a,T)  and  R^(a,T)  are  tp(a,T)  =  exp[yp(0,T)],  where 
yp(a,T)  =  a(a,T)  Wp  +  b(a,T),  and  R^(a,T)  =  exp{-exp[(s,n  t  -  b(a,T))/a(a,T)]} . 
Estimated  standard  errors  of  these  ML  estimates  may  be  computed  from  results 
given  in  Appendix  B. 


Plots  of  selected  estimated  lifetime  quantiles  are  presented  in 
Figures  6  through  13.  In  Figure  6,  the  ML  estimates  tp(a,T)  are  displayed  for 
the  quantile  probability  p  =  .5  (representing  the  median),  the  temperatures 
65°C,  45°C,  25°C,  and  the  stress  range  200  to  450  ksi.  A  similar  plot  for  the 
quantile  probability  p  =  10"^  is  presented  in  Figure  7.  In  Figures  8  and  9, 


the  estimates  tp(a,T)  are  shown  for  the  quantile  probabilities  p  =  .5  and 
-3 

p  =  10  ,  respectively,  seven  selected  stress  levels,  and  the  temperature 
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range  20°  to  70°C.  In  Figures  10  and  11,  the  temperature  is  held  fixed,  at 
25°  and  45°C,  respectively,  and  the  estimates  tp(a, t)  are  displayed  for 
selected  quantiles  and  a  stress  range.  Finally,  in  Figures  12  and  13,  the 
stress  level  is  held  fixed,  at  100  ksi  and  200  ksi,  respectively,  and  the 
estimates  tp(a,  t)  are  plotted  for  selected  quantiles  and  a  temperature 
range . 


4.  Discussion 

The  sources  of  uncertainty  (a)  through  (d)  given  in  part  A,  with  the 
obvious  extension  in  (c)  from  stress  a  to  (stress,  temperature)  =  (a, t),  are 
valid  in  the  present  context  as  well.  In  addition,  there  are  uncertainties 
involving  (e)  the  unquantified  (assumed  negligible)  UV  affect  at  the  four  room 
temperature  stress  levels;  (f)  possible  (assumed  nonexistent)  physical 
differences  between  the  strands  used  for  room  temperature  testing  and  the 
strands  used  for  elevated  temperature  testing;  and  (g)  the  unquantified  and 
unmodeled  effect  of  a  higher  relative  humidity  at  room  temperature  than  at 
elevated  temperatures. 

C.  NASA  Kevlar  49/epoxy  spherical  pressure  vessel  data  analysis:  lifetime 
versus  pressure  (room  temperature  environment  without  UV  light). 

1.  Materials 

The  Kevlar  49/epoxy  vessels  used  in  room  temperature  strength  and 
lifetime  tests  are  described  in  Table  5.  Added  detail  and  descriptions  of  the 
testing  equipment  may  be  found  in  References  3,  5,  7,  and  8. 
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Table  5.  Summary  of  material  characteristics.  (Room  temperature  vessels) 

material  system:  Kevlar  49  (380  denier,  approximately  267  filaments)  wound 

in  Dow  DER  332/Jefferson  Jeff amine  T-403  with  weight  ratio 
100/44.7 

vessel  geometry:  112  mm  interior  diameter  spherical  aluminum  liner,  1  mm 

thick,  overwrapped  with  1.1  mm  composite  wall,  wound  in  a 
delta  axisymmetric  pattern. 

2.  Data 

The  NASA  pressure  vessels  were  separated  into  five  stress 
environments  for  stress-rupture  testing  at  room  temperature  without  light. 
Initial  burst  tests  resulted  in  an  estimated  static  strength  of  5.0  ksi.  The 
complete  set  of  burst  loads  is  presented  in  Appendix  C.  Each  vessel  was  wound 
from  one  of  eight  different  Kevlar  spools  of  yarn.  It  was  subsequently 
discovered  (Reference  3)  that  one  of  the  spools  (spool  #7)  had  physical 
properties  distinctively  different  from  the  others  (in  terms  of  filament  number 
and  diameter,  sodium  hydroxide  content,  etc.).  It  was  appropriate,  then,  to 
remove  from  consideration  all  data  pertaining  to  specimens  from  this  spool. 

The  statistical  analysis  here  is  based  on  the  exact,  grouped,  and  censored 
failure  times  of  the  139  vessels  wound  from  the  remaining  seven  spools, 
numbered  1  through  7  for  convenience.  (Hence  spool  #7  in  this  report  refers 
to  the  original  spool  #8.)  The  stress-rupture  data  are  summarized  in  Table  6. 
Spool  identity  is  included  in  the  classification  due,  as  noted  in  the  next 
section,  to  a  significant  spool  effect. 
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3.  Statistical  Analysis  of  the  Data 

Again  the  basic  Weibull  lifetime  distribution  asumption  is  supported 
by  Weibull  probability  paper  plots.  An  initial  investigation  of  the  data 
ignored  possible  spool  effects  and  treated  all  times  at  a  given  stess  level  as 
a  sample  from  a  single  Weibull  distribution.  From  the  marginal  ML  estimates,  (a, 
b)  =  (1.75,  5.51),  (1.14,  6.73),  (.952,  9.01),  and  (.935,  11.2)  for  the  respective 
pressure  levels  4.28,  3.97,  3.68,  and  3.38  ksi,  the  dependencies  a (a)  =  l/a(a) 
and  b(a)  =  £n  3(a)  appeared  to  be  modeled  adequately  by 

(5)  a(a)  =  6^  +  ^2^^  b(a)  =  0^  +  +  Q^(P 

Subsequently,  in  response  to  work  by  Gerstle  (Reference  3)  who  drew  attention 
to  spool-to-spool  variability,  it  was  decided  to  incorporate  spool  effects 
into  the  dependency  models.  This  was  done  by  allowing  the  transformed  scale 
parameter  b  to  depend  upon  spool  identity.  The  model  ultimately  selected  was 

(6)  a(a)  =  +  S^e^*^  and  b(a)  =  +  0^6^  +  6^6^ 


+  SgSg  +  a  , 

where  the  <5^^  are  indicators  of  spool  identity: 

^  1  if  the  vessel  is  wound  from  spool  i 

=  Jo  if  the  vessel  is  wound  from  a  spool  other  than  spool  i  or  spool  7 

[-1  if  the  vessel  is  wound  from  spool  7, 

i  =  1,  ...,  6.  The  effect  of  spool  i  for  1  _<  i  _<6  is  effect  of 

spool  7  is  -(0^  +  ...  +  0q),  and  the  sum  of  all  spool  effects  is  zero.  Thus, 
spool-to-spool  variation  is  parametrized  in  terms  of  relative  additive  factors 
in  transformed  scale.  If  0^  equaled  0.60,  say,  then  vessels  wound  from 
spool  3  would  have  lifetimes  with  scale  parameter  (characteristic  life) 
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stress. 

Table 

{%  UTS) 

6.  Room  temperature  vessel  data 

spool#  #vessels  #exact 

summary. 

#grouped 

#censored 

ksi 

4.28 

(85.6) 

1 

4 

failure 

times 

4 

failure 

times 

0 

failure 

times 

0 

4.28 

2 

8 

8 

0 

0 

4. 28 

3 

4 

4 

0 

0 

4.28 

4 

5 

5 

0 

0 

4. 28 

5 

4 

4 

0 

0 

4.28 

6 

3 

3 

0 

0 

4.28 

7 

3 

3 

0 

0 

3.97 

(79.4) 

1 

4 

4 

0 

0 

3.97 

2 

7 

7 

0 

0 

3.97 

3 

4 

4 

0 

0 

3.97 

4 

4 

4 

0 

0 

3.97 

5 

1 

1 

0 

0 

3.97 

6 

2 

2 

0 

0 

3.97 

7 

2 

2 

0 

0 

3.68 

(73.6) 

1 

3 

3 

0 

0 

3.68 

2 

6 

6 

0 

0 

3.68 

3 

2 

2 

0 

0 

3.68 

4 

2 

2 

0 

0 

3.68 

5 

1 

1 

0 

0 

3.68 

6 

3 

3 

0 

0 

3.68 

7 

6 

5 

1 

0 

3.38 

(67.6) 

1 

4 

0 

0 

4 

3.38 

2 

1 

1 

0 

0 

3.38 

3 

2 

1 

0 

1 

3.38 

4 

5 

0 

0 

5 

3.38 

5 

3 

3 

0 

0 

3.38 

6 

3 

3 

0 

0 

3.38 

7 

3 

0 

0 

3 

2.49 

(49.8) 

1 

5 

0 

0 

5 

2.49 

2 

6 

0 

0 

6 

2.49 

3 

3 

0 

0 

3 

2.49 

4 

7 

0 

0 

7 

2.49 

5 

8 

0 

0 

8 

2.49 

6 

7 

0 

0 

7 

2.49 

7 

4 

0 

0 

4 

,0.60  , 

1.82  times  as 

great  as 

the  geometric 

mean  characteristic 

life  among 

the  seven  spools.  The  transformed  shape  parameter  a,  on  the  other  hand,  is 


assumed  to  depend  only  on  stress  and  not  on  spool  identity. 


The  ML  fit  obtained  by  assuming  the  model  (6),  and  applying  to  the  failure 
time  data  the  basic  methodology  of  analyses  A  and  B,  showed  a  substantial  and 
dramatic  improvement,  as  measured  by  maximized  log  likelihood,  over  the 
analogous  ML  fit  based  upon  the  model  (5)  which  does  not  consider  any  spool 


effects.  The  ML  estimates  of  the  coefficients  (e^^,  ...,  are  the 
following:  (0^,  02,  ©3,  0^,  ©5,  ©g,  0^,  0g,  0^g)  =  (3.20037  x  10"^, 

1.51037  X  10"^,  1.09470,  -8.63305  x  10'^,  -1.58238,  1.49728,  -1.48675  x  10"^, 
-4.28710  X  10”^,  3.15136  x  10^,  -6.18686).  The  ML  estimated  effect  of  spool  7 
is  -(©3  +  ...  +  0q)  =  4.31090  X  10"^.  The  ML  estimates  of  the  functions  a,  b, 

^  A  ^  Orr  ^  .A  ^  /S  A  A 

a  and  ©  are  a(a)  =  ©^  +  ©2e  ,  b(a)  =  ©j&j^  +  0^^  +  ©^&j  +  0^6^  +  0^^  +  0g6g  + 

A  A  ^  A  ^ 

©9  +  =  l/a(o),  and  B(a)  =  exp[b(a)].  In  addition,  the  ML  estimates 


of  the  lifetime  parameters  t  (a)  and  R^(a)  are  t  (a)  =  exp[y  (a)],  where  y^Ccf) 

pup  P  P 

=  a  (a)  Wp  +  b(a),  and  R|.(ct)  =  expi-exp[(Sn  t  -  b(a))/a(a)] }.  For  inferences 
pertinent  to  the  "average”  or  "typical"  spool,  i.e.,  the  conceptual  spool 
corresponding  to  a  zero-valued  spool  effect,  ML  estimates  of  the  above 
quantities  are  obtained  by  substituting  63  =  . . .  =  6g  =  0  into  the  expressions 

A  A  A  A. 

for  b  and  b.  In  this  case  b(a)  =  9^  +  9^^^  a  estimates  b(a)  =  9^  +  a. 
Estimated  standard  errors  of  all  ML  estimates  are  computable  from  results 
presented  in  Appendix  C. 

The  maximum  likelihood  estimates  of  selected  lifetime  parameters  as 
functions  of  spool  and  stress  are  displayed  in  the  accompanying  Figures  14 
through  17.  In  Figure  14  and  15,  ML  estimates  of  quantiles  and  reliabilities, 
respectively,  are  plotted  for  the  "average"  spool  case.  In  Figures  16  and  17, 
the  magni’tude  of  spool-to-spool  variation  is  illustrated  in  plots  of  the 
estimated  10“^-quantiles  and  5-year  reliabilities,  respectively,  for  the 
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"best,"  "average,"  and  "worst"  spools.  The  spool  which  performed  "best"  in  the 
testing  was  spool  4,  whose  estimated  spool  effect,  0g  =  1.49728,  was  the 
largest  among  the  seven  spools.  The  "worst"  spool,  corresponding  to  the 
smallest  estimated  spool  effect  (e^  =  -1.58238),  was  spool  3.  Thus  the  "best" 
estimates  are  based  0063=64  =  65=67=63  =  0,  6g  =  1;  the  "worst" 
estimates  0063  =  64=65  =  67  =  63  =  0,  65  =  1;  and  the  "average"  estimates 
on  63  =  64  =  65  =  65  =  57  =  63  =  0. 

4.  Discussion 

The  sources  of  uncertainty  (a)  through  (d)  inroduced  in  part  A  for 
strands  are  also  pertinent  in  the  vessel  context.  An  added  concern  here 
involves  the  relevance  of  the  particular  spools  used  in  the  experiment. 

Ideally  these  seven  spools  constitute  a  random  sample  of  the  conceptual 
population  of  all  spools  producible  and  suitable  for  construction  of  vessels. 
Estimates  based  on  the  "average"  spool  effect  (obtained  above  by  setting  each 
6 ^  =  0,  i  =  3,  . . . ,  8)  then  pertain  to  the  population  average  spool,  a 
conceptual  spool  of  interest.  In  addition,  a  measure  of  the  dispersion  of 
spool  effects  within  the  conceptual  spool  population  is  available,  namely  the 
sample  standard  deviation,  1.086,  of  the  seven  ML  estimated  spool  effects. 
Unfortunately,  however,  the  seven  spools  appear  to  have  been  selected  by 
convenience  rather  than  by  a  random  scheme.  Consequently,  inferences  cannot 
properly  be  extended  to  spools  beyond  the  seven  tested,  and  the  "average" 
spool  results  pertain  merely  to  a  conceptual  spool  which  typifies  the  seven. 


D.  Power  law  fits 


In  the  absence  of  sources  of  chemical  degradation,  such  as  UV  light  or 
humidity  or  elevated  temperatures,  it  has  been  argued  (see  Reference  1)  that 
the  power  law  model  for  lifetimes  holds.  Thus  for  Kevlar  A9/epoxy  strands  at 
room  temperature  in  darkness  the  model 

a(o)  =  and  b(a)  =  02  +  6^  itn  0 

would  be  appropriate.  Based  on  the  failure  time  data  corresponding  to  the 
upper  four  stress  levels  from  the  table  of  part  A  (for  which  the  UV  effect  may 
be  assumed  negligible  or  slight),  the  ML  methodology  yields  the  estimates 

6^^  =  1.04163  ,  02  =  2.78954  x  10^,  and  0^  =  -4.55850  x  10^  . 

Figure  18  displays  corresponding  ML  estimates  of  selected  quantiles.  The 
results  are  theoretically  appropriate  for  strands  under  conditions  of  no 
chemical  degradation. 

The  power  law  fit  applied  to  the  NASA  pressure  vessel  data  (motivated  by 
experimental  conditions  of  darkness  and  room  temperature),  does  not  yield  ML 
results  competitive  in  terms  of  maximized  log  likelihood  to  the  fit  obtained 
using  model  (6)  in  part  C.  (In  fairness,  the  aforementioned  theoretical  work 
applies  to  strands  but  not  to  structures.)  Nonetheless,  the  power  law  model 

a(a)  =  0j^  and  b(o)  =  02*^1  ®4'^3  ®5^4  ®6%  ®8  ^  is 

estimated  by  (0^^,  ®6»  ®7'  ®8’  ®9^  ^  (7.81416  x  10'^,  1.17832, 

-9.44661  X  10"^,  -1.65797,  1.68801,  -3.21649  x  10"^,  -5.11371  x  10"^,  3.92902  x 

10^,  -2.34652  X  10^).  The  ML  estimated  spool  7  effect  is  therefore  -(©2  +  ...  + 
-1 

0^)  =  5.69321  X  10'  .  In  Figure  19,  the  power  law  ML  estimates  of  selected 
quantiles  are  plotted  in  terms  of  stress  for  the  "average”  spool. 
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E.  Appendix  A 

Table  7  gives  the  ordered  breaking  loads  for  53  room  temperature  strand 
specimens. 

Tables  8  through  14  give  the  observed  lifetimes  for  room  temperature 
strands  at  seven  stress  levels. 

Formulas  presented  here  allow  straightforward  computation  of  ML  estimates 
and  their  estimated  standard  errors  (i.e.,  square  roots  of  estimated  variances) 
for  quantities  introduced  in  part  A. 

From  equations  (2)  and  (3),  it  follows  that  y_(a)  can  be  expressed  as 
a  third  degree  polynomial  in  a,  namely 

^  <*•  A  A  AQAA 

yp(a)  =  (e^  Wp  +  64)  +  (02  Wp  +  65)  0  +  (63  Wp  +  65)  a  +  67  o  . 

A 

The  graphed  versions  of  Figures  1,  3,  and  4  involve  log  tp(CT)  =  (log  e) 
p(a),  where  log  e  =  0.4342944819. 


Table  7.  Ordered  breakijig  loads,  in 
49/epoxy  strand  specimens. 

Breaking  load,  lb. 

20.5 

20.9 

21.1 

21.4 

21.6 
21.8 
22.0 
22.2 

22.5 
22.7 

22.9 
23.1 


pounds,  of  room  temperature  Kevlar 

Quantity 

1 

1 

1 

5 

6 

4 

13 

5 
8 
5 
3 

_i 

53 


-24- 


TABLE  8. 

Ranks  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Strands 


Stress:  455.2  ksl 
Temperature:  25  deg  C. 
Number  of  specimens:  101 


1 

9 

27 

0 

.24 

53 

0.83 

79 

1 

.58 

2 

9 

28 

0 

.29 

54 

0.85 

80 

1 

.53 

3 

0 

.02 

29 

0 

.34 

55 

0.90 

81 

1 

.54 

4 

0 

.02 

30 

0 

.35 

56 

0.92 

88 

1 

.54 

5 

0 

.02 

3.1 

0 

.36 

57 

0.95 

83 

1 

.55 

6 

0 

.03 

32 

0 

.38 

58 

0.99 

84 

1 

.58 

7 

0 

.03 

33 

0 

.40 

59 

1  .00 

85 

1 

.60 

8 

0 

.04 

34 

0 

.42 

60 

1  .01 

86 

1 

.63 

9 

0 

.05 

35 

0 

.43 

61 

1 .02 

87 

1 

.64 

10 

0 

.06 

36 

0 

.52 

62 

1 .03 

88 

1 

.80 

1  1 

0 

.07 

37 

0. 

.54 

63 

1 .05 

89 

1  , 

.80 

12 

0 

.07 

38 

0, 

.56 

64 

1  .  10 

90 

1  , 

.81 

13 

0 

.08 

39 

0, 

.60 

65 

1.10 

91 

2, 

.08 

14 

0 

.09 

40 

0. 

.60 

66 

1.11 

98 

2, 

.05 

15 

0  , 

.09 

41 

0. 

.63 

67 

1  .  15 

93 

2, 

.  14 

16 

0, 

.  10 

42 

0. 

.65 

68 

1  .  18 

94 

2. 

.  17 

17 

0, 

.  10 

43 

0, 

.67 

69 

1  .20 

95 

.33 

18 

0. 

.  1  1 

44 

0, 

.68 

70 

1.29 

96 

3. 

.03 

19 

0. 

.  1  1 

45 

0. 

,72 

71 

1  .31 

97 

3. 

,  03 

20 

0. 

,  12 

46 

0. 

,72 

72 

1  .33 

98 

3. 

,34 

21 

0, 

,  13 

47 

0. 

.72 

73 

1 .34 

99 

4. 

,80 

22 

0, 

,  18 

48 

0. 

73 

74 

1  .40 

100 

4, 

69 

23 

0. 

,  19 

49 

0. 

79 

75 

1  .43 

101 

7. 

89 

24 

0. 

20 

50 

0, 

79 

76 

1  .45 

25 

0. 

23 

51 

0. 

80 

77 

1  .50 

26 

0. 

24 

52 

0. 

80 

78 

1 .51 

Grouped 

time 

in  the 

i  n 

1 1  e  r  V  a  1 

(0. 

00,0.01  )  . 

25 


TABLE  9. 

RanUs  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Strands 


Stress:  440 .0  Us\ 
Temperature:  25  deg  C. 
Number  of  specimens:  100 


1 

0.03 

26 

1 .35 

51 

o 

00 

ru 

76 

5.92 

2 

rr 

o 

o 

27 

1 .35 

52 

2.89 

77 

6.19 

3 

o 

o 

28 

1 .36 

53 

2.91 

78 

6.33 

4 

0.08 

29 

1.39 

54 

2.91 

79 

6.59 

5 

0.08 

30 

1 .40 

55 

3. 10 

80 

ru 

UD 

CD 

6 

0.09 

31 

1  .48 

56 

3. 1  1 

81 

CD 

CD 

7 

0.  15 

32 

1  .55 

57 

3.42 

82 

7. 16 

8 

0.  15 

33 

1 .56 

58 

3.47 

83 

7.29 

9 

0.20 

34 

1  .57 

59 

3.50 

84 

7.36 

10 

0.28 

35 

1  .59 

60 

3.64 

85 

7.38 

1  1 

0.28 

36 

1 .80 

61 

3.67 

86 

7.42 

12 

0.30 

37 

1 .85 

62 

3.68 

87 

8.09 

13 

0.38 

38 

1  .86 

63 

3,75 

88 

8.88 

14 

0.41 

39 

1  .92 

64 

4.39 

89 

9.05 

15 

0.50 

40 

2.00 

65 

4.50 

90 

9.25 

16 

0.51 

41 

2.08 

66 

4.53 

91 

9.26 

17 

0.57 

42 

2.27 

67 

4.53 

92 

9.47 

18 

0.70 

43 

2.38 

68 

4.65 

93 

10.50 

19 

0.72 

44 

2.46 

69 

4.65 

94 

10.57 

20 

ru 

00 

o 

45 

2.49 

70 

4.69 

95 

11.25 

21 

0.84 

46 

2.61 

71 

4.70 

96 

13.49 

22 

1 .05 

47 

2.61 

72 

4.75 

97 

13.78 

23 

1.13 

48 

2.62 

73 

4.84 

98 

16.  15 

24 

1  .  19 

'  49 

2.74 

74 

5.01 

99 

16.59 

25 

1 .32 

50 

ru 

CD 

75 

5.71 

100 

18.  16 

26 


TABLE  10, 

Ranks  and  Breaking  Times  (hr) 
Stress:  HB4 .8  k  s I 


Temperature:  85 

deg 

C. 

Number 

of  specimens: 

;  103 

1 

0.25 

27 

7.32 

2 

0.31 

28 

7.45 

3 

0.44 

29 

CD 

ru 

4 

0.45 

30 

7.99 

5 

0.57 

31 

CD 

ru 

ru 

6 

0.65 

32 

8.35 

7 

0.73 

33 

8.38 

8 

0.89 

34 

8.39 

9 

1 .06 

35 

8.45 

10 

1  .22 

36 

8.53 

1  1 

1 .37 

37 

8.55 

12 

1.83 

38 

8.64 

13 

1.96 

39 

8.68 

14 

2.15 

40 

00 

CD 

ru 

15 

2.40 

41 

8.93 

16 

2.51 

42 

9.45 

17 

2.77 

43 

9.57 

18 

4.05 

44 

9.80 

19 

4.07 

45 

9.83 

20 

4.34 

46 

10.60 

21 

4.98 

47 

10.82 

22 

5.86 

48 

10.83 

23 

5.90 

49 

11.03 

24 

6.18 

50 

11.12 

25 

6.30 

51 

11.13 

26 

7.14 

52 

12.52 

for  Kevlar  49/Epoxy  Strands 


53 

12 

.94 

79 

34 

.51 

54 

16 

.42 

80 

34 

.74 

55 

16 

.47 

81 

35 

.63 

56 

17 

00 

CD 

82 

36 

.84 

57 

17 

.79 

83 

36 

.97 

58 

18 

.45 

84 

39 

.  17 

59 

18 

.95 

85 

41 

.54 

60 

19 

.07 

86 

42 

.32 

61 

19 

.40 

87 

42 

.53 

62 

19 

.62 

88 

43 

.61 

63 

19 

.86 

89 

48, 

.54 

64 

20, 

.76 

90 

49, 

.02 

65 

21  . 

.38 

91 

55. 

.20 

66 

23. 

.03 

92 

55. 

CD 

CD 

67 

23. 

.  10 

93 

61  , 

.37 

68 

23. 

83 

94 

63. 

,  17 

69 

24. 

46 

95 

66. 

,48 

70 

24. 

81 

96 

67. 

06 

71 

25. 

15 

97 

74. 

01 

72 

25. 

18 

98 

74. 

61 

73 

26. 

96 

99 

76. 

46 

74 

27. 

53 

100 

84. 

26 

75 

27. 

86 

101 

89. 

87 

76 

29. 

89 

102 

97. 

37 

77 

32. 

55 

103 

1  19. 

09 

78 

33. 

95 

27 


TABLE  I  1  . 

Ranks  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Strands 

Stress:  404 .6  k  s  1 
Temperature:  25  deg  C. 


Number 

0  f 

spec 

i me  ns  : 

100 

1 

1 

.8 

26 

84, 

.2 

51 

152 

.2 

76 

285 

.9 

2 

3 

.  1 

27 

87, 

.  1 

52 

152, 

.8 

77 

292 

.6 

3 

4 

.2 

28 

87, 

.3 

53 

157, 

.7 

78 

295, 

.  1 

4 

6 

.0 

29 

93. 

.2 

54 

160. 

,0 

79 

301  , 

.  1 

5 

7 

.5 

30 

103, 

.4 

55 

163. 

.6 

80 

304, 

.3 

6 

8 

.2 

31 

104. 

.6 

56 

166. 

,9 

81 

316, 

.8 

7 

8 

.5 

32 

105. 

.5 

57 

170, 

,5 

82 

329, 

.8 

8 

10 

.3 

33 

108. 

,8 

58 

174. 

,9 

83 

334, 

,  1 

9 

10 

.6 

34 

1  12, 

.6 

59 

177. 

.7 

84 

346, 

,2 

10 

24 

.2 

35 

1 16, 

.8 

60 

179. 

,2 

85 

351  , 

,2 

1  1 

29 

.6 

36 

1 18. 

.0  ■ 

61 

183. 

,6 

86 

353, 

.3 

12 

31 

.7 

37 

122. 

.0 

62 

183. 

,8 

87 

369. 

.3 

13 

41 

.9 

38 

123. 

.5 

63 

194, 

.3 

88 

372, 

,3 

14 

44 

.  1 

39 

124, 

.4 

64 

195. 

,  1 

89 

381  , 

.3 

15 

49 

.5 

40 

125, 

.4 

65 

195, 

.3 

90 

393. 

.5 

16 

50 

.  1 

41 

129. 

.5 

66 

202. 

.6 

91 

451  . 

.3 

17 

59 

.7 

42 

130, 

.4 

67 

220  . 

.2 

92 

461  . 

,5 

18 

61 

.7 

43 

131  . 

.6 

68 

221  . 

,3 

93 

574, 

,2 

19 

64 

.4 

44 

132. 

.8 

69 

227, 

.2 

94 

653. 

,3 

20 

69 

.7 

45 

133, 

.8 

70 

251  , 

.0 

95 

663, 

.0 

21 

70 

.0 

46 

137, 

.0 

71 

266. 

.5 

96 

669. 

,8 

22 

77 

.8 

47 

140. 

,2 

72 

267, 

.9 

97 

739. 

,7 

23 

80 

.5 

48 

140. 

.9 

73 

269. 

,2 

98 

759. 

.6 

24 

82 

,3 

49 

148. 

.5 

74 

270, 

,4 

99 

894, 

.7 

25 

83 

.5 

50 

149. 

.2 

75 

272, 

,5 

100 

974, 

,9 

28 


TABLE  12. 


Ranks 

and  Break! 

ng 

Times  ( hr ) 

for 

Kevlar 

49/Ep  oxy  Strands 

St  r  e  s.s 

:  35H.0  ksi 

Temperature:  25 

deg 

C. 

Number 

of  specimens 

:  49 

1 

1051. 

14 

5817. 

27 

9711. 

40 

12044. 

2 

1337. 

15 

5905. 

28 

9806. 

41 

13520. 

3 

1389. 

16 

5956. 

29 

10205. 

42 

13670. 

4 

1981  . 

17 

6068. 

30 

10396. 

43 

14110. 

5 

19H2. 

18 

6121  . 

31 

10861 . 

44 

14496. 

6 

2328. 

19 

6473. 

32 

11026. 

45 

15395. 

7 

3629. 

20 

7501  . 

33 

1 1214. 

46 

16179. 

8 

4006. 

21 

7886. 

34 

1 1362. 

47 

17092. 

9 

4012. 

22 

8108. 

35 

1 1604. 

48 

9 

10 

4063. 

23 

8546. 

36 

1 1608. 

49 

9 

1 1 

4921  . 

24 

8666. 

37 

1  1745. 

12 

5445. 

25 

8831  . 

38 

1 1762. 

13 

5620. 

26 

9106. 

39 

1 1895. 

n 


g  Grouped  time 


the  interval 


(17408. , 17576. )  . 


TABLE  13. 


Ranks  and  Break! 

ing  Times  (hr) 

for 

Ke  v  1  a 

Stress:  303 .5  h  s i 

Temperature:  25 

deg  C . 

Number  of  specimens:  50 

1  1387a. 

14  89838. 

87 

39840 

a  i8oa4. 

15  31884. 

88 

39576 

3  19008. 

16  31758. 

89 

39744 

4  ai960. 

17  38838. 

30 

39744 , 

5  aa87a. 

18  38976. 

31 

41598, 

6  85008. 

19  35544. 

38 

41760, 

7  85848. 

80  35760. 

33 

41760. 

8  87816. 

81  35988. 

34 

48600. 

9  87744. 

88  36588. 

35 

48600. 

10  87840. 

83  36780. 

36 

48960 . 

11  88518. 

84  38598. 

37 

43176. 

18  88896. 

25  38592. 

38 

44664 . 

13  89838. 

26  39078. 

39 

49176. 

H9/Ep  0 X  y  SI r  a n 


40  494H0 

41  5119a 

42  gl 

43  gl 

44  ga 

45  ga 

46  ga 

47  ga 

48  ga 

49  ga 

50  ga 


gl  Grouped  time  In  the  interval  ( 49080 66408 . 
ga  Grouped  time  in  the  interval  ( 54744 7a07a .) . 


r 


TABLE  14. 

Ranks  and  Breaking  Timfes  (hr)  for  Kevlar 

Stress;  252 .9  k  s i 
Temperature:  25  deg  C. 

Number  of  specimens:  50 


1  3134H.  2  32376.  3  58056. 


Censored  Times  and  Quantities 

10000.  1 
14376.  1 

66408.  9 

72072.  35 


49/Epoxy  Strands 


4  66024. 


31 


From  large  sample  statistical  theory,  the  joint  probability  distribution 
of  the  ML  estimates  (0i,  67)  is  approximately  7-dimensional  multivariate 
normal  with  mean  vector  (0p  .,.,07)  and  covariance  matrix  Z  estimated  by 

^  ^  ***’ 


4.7576  X  10 
-2.8951  X  10 
4.2251  X  10 
2.1477  X  10*^ 
-1.9355  X  10 
5.7394  X  10 


-1 


-6 


-2 


-5 


-2.8951  X  10 


-3 


4.2251  X  10 


-6 


2.1477  X  10°  -1.9355  x  10 


-2 


5.7394  X  10" 


-5.5780  X  10 


1.7743  X  10'^  -2.6046  x  10"°  -1.2800  x  10"^  1.1607  x  10"^  -3.4605  x  10"^  3.3785  x  10 


-2.6046  X  10"°  3.8449  x  10"^^  1.8434  x  10"°  -1.6787  x  10"^ 


5.0227  X  10 


■10 


-4.9186  x  10 


-8 

-10 

■13 


-1.2800  X  10 


-2 


1.8434  X  10 


4.2525  X  10^  -3.7848  x  10"^  1.1021  x  10"°  -1.0494  x  10 


-6 


1.1607  X  10"^  -1.6787  x  10"^  -3.7848  x  10"^  3.3820  x  10"° 


-9.8833  X  10 


-6 

^-8 


9.4416  X  10" 


-11 


,-8 


-3.4605  X  10"^  5.0227  x  10“^°  1.1021  x  10"°  -9.8833  x  10"°  2.8980  x  10  “  -2.7770  x  10 

X  10"’°  -4.9186  X  10"’°  -1.0494  X  10"°  9.4416x10"°  -2.7770  x  10"’’  -2.6685  x  10"’"’ 


The  corresponding  estimates  of  the  variances  of  the  estimators  a(a)  and  b(a) 

2 

are  expressible  as  polynomials  in  a  (where  (c-j,  ...  ,  Cy)  -  (1,  a»  a  >  li 

a,  a^,.a^)) : 

^  .  3  3 


Var(a(a))  =  c.  c.  s..  =  +  2 


a  (S22  +  2  5^3)0 


+  2  $23  a  +  ® 


and 


Ab{a))  =  J  I  =  ^44  *  ^  ^45  ^*55  ^  ^  ^6^'^' 


i=4  j=4 
M7 


+  2(s^7  +  ^56^^  *  ^^66  ^  ^57^°  ’’’  ^  ^67 


’66 


77 
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A  ^ 

Similarly,  the  estimated  covariance  between  a(a)  and  b(a)  is  expressible  as 

3  7 

b(o))  ^  C,CjS,j  >  s,^  +  (s,5  *  t  (s,5  t  +  53^)0  2 

*^17  *  ^26  *  *35*''  *  ^*27  *  *37  " 

A  ^ 

Thus,  estimates  of  the  variances  of  y  (a)  and  R^(a)  are 

♦  p  U 

^(yp(a))  =  Wp^  v^(a(a))  +  V^(b(a))  +  2  Wp  C^(a(a),  b(a)),  and 

Var(L(a))  =  |J —  exp[(jin  t  -  b(a))/a(a)]  R^(0) 

^  |a(a)  ^ 

{[(jin  t  -  b(a))/a(a)]^  '^(a(a))  +  V^^(b(a)) 

+  2[(jin  t  -  b(a))/a(a)]  {a(a),  b(CT))}  . 

A  confidence  interval  for  the  quantile  yp(a)  with  confidence  coefficient 

approximately  1-a  is  therefore  y  (a)  +  z,  ,oC'^(y  where  z 

p  —  i-a/z  p  Y 

denotes  the  y-quantile  of  the  standard  normal  distribution  (e.g.,  z  95  =  1.645). 
The  corresponding  confidence  interval  for  the  quantile  tp(a)  is  accordingly 

Cexp{yp(a)-z^^/2ttf^^(yp(«^))^^^^}»  exp{yp(a)  +  Zi^/2l^^^yp(a))]^^^}]. 

Similarly,  a  lower  confidence  bound  for  Rj.(cf)  with  confidence  coefficient 
approximately  1-a  is  R^(a)  -  z^_^[V^(R^(a))]^^^. 

F.  Appendix  B 

Results  analogous  to  those  of  Appendix  A  are  presented  here  to  provide 
assistance  in  computations  of  quantities  introduced  in  part  B. 
fi  Strength  tests  were  performed  on  13  strand  specimens  from  the  lot  used  in 

subsequent  elevated  temperature  stress-rupture  tests.  For  the  13  strands. 
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breaking  loads  were  measured  at  room  temperature.  The  Individual  values  have 
been  lost.  However,  their  arithmetic  average,  22.66  lb,  has  survived. 

Tables  15  through  27  give  the  lifetime  data  for  the  13  environments  at 
elevated  temperatures.  (The  lifetime  data  for  the  four  environments  at  room 
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TABLE  15. 

Ranks  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Strands 


Stress:  436.7  ksi 
Temperature:  180  deg  C. 
Number  of  specimens:  48 


1 

0.0068 

13 

0.0162 

a 

0.007P 

14 

0.0175 

3 

0.0080 

15 

0.0177 

H 

0.0095 

16 

0.0182 

5 

0.0118 

17 

0.0185 

6 

0.0122 

18 

0.0190 

7 

0.0130 

19 

0.0217 

8 

0.0132 

20 

0 . 0233 

9 

0.0150 

21 

0.0233 

10 

0.0152 

22 

•0 . 0235 

1  1 

0.0152 

23 

0 . 0242 

IE 

0.0162 

24 

0 . 0245 

25 

0.0247 

37 

0.0425 

26 

0.0262 

38 

0.0478 

27 

0.0290 

39 

0.0485 

28 

0.0308 

40 

0 . 0547 

29 

0.0317 

41 

0 . 0665 

30 

0.0332 

42 

0.0700 

31 

0.0333 

43 

0.0747 

32 

0.0342 

44 

0 . 0752 

33 

0.0350 

45 

0.0763 

34 

0.0383 

46 

0.0815 

35 

0.0403 

47 

0.0830 

36 

0.0413 

48 

0. 1247 

35 


TABLE  16. 

Ranks  and  Breaking  Times  (hr)  for  Kevlar  H9/Epoxy  Strands 


Stress :  436.7  ks i 
Temperature:  110  deg  C. 
Number  of  specimens:  39 


1 

0.0148 

1  1 

0.0593 

2 

0.0248 

12 

0.0595 

3 

0.0263 

13 

0.0605 

4 

0.0301 

14 

0.0618 

5 

0.0323 

15 

0.0713 

6 

0.0461 

16 

0.0728 

7 

0.0530 

17 

0.0740 

8 

0.0573 

18 

0.0770 

9 

0.0580 

19 

0.0775 

10 

0.0593 

20 

0.0830 

21 

0.0885 

31 

0.1395 

22 

0.0900 

32 

0.  1410 

23 

0.0910 

33 

0. 1448 

24 

0.0928 

34 

0. 1516 

25 

0.0933 

35 

0.1663 

26 

0.0963 

36 

0. 1727 

27 

0.1113 

37 

ru 

CD 

o 

28 

0. 1240 

38 

0. 1783 

29 

0.1330 

39 

0. 1783 

30 

0.1347 

36 


TABLE  17, 

Panics  and  Braakihg  Times  (hr)  for  Kevlar  49/Epoxy  Strands 

Stress  :  436.7  ks i 
Temperature:  100  deg  C. 

Number  of  specimens:  46 


1 

0.0407 

13 

0.0890 

25 

0.1517 

37 

0.2080 

2 

0.0437 

14  ■ 

0 . 0977 

26 

0.1567 

38 

0.2249 

3 

0.04B3 

15 

0.1063 

27 

0.1612 

39 

0 . 2262 

4 

0.0480 

16 

0 . 1 087 

28 

0.1613 

40 

0.2432 

5 

0.0526 

17 

0.1119 

29 

0.1662 

41 

0.2478 

6 

0.0532 

18 

0.1223 

30 

0.1750 

42 

0 . 2497 

7 

0.0545 

19 

0.1233 

31 

0.1876 

43 

0 . 2537 

8 

0.0587 

20 

0.1347 

32 

0.1912 

44 

0.2873 

9 

0.0697 

21 

0 . 1 362 

33 

0.1932 

45 

0.2917 

10 

0.0718 

22 

0.1432 

34 

0.1950 

46 

0 . 3762 

1 1 

0.0770 

23 

0.1463 

35 

0.1963 

12 

0.0837 

24 

0. 1500 

36 

0.2032 

37 


TABLE  18. 

Ranks  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Strands 


Stress:  410.7  ksl 
Temperature:  120  deg  C. 
Number  of  specimens:  46 


1 

O.OIOE 

13 

0 . 0997 

2 

o.oaoo 

14 

0.1078 

3 

0.0200 

15 

0. 1082 

4 

0.0612 

16 

0.1113 

5 

0.0650 

17 

0.1170 

6 

0.0738 

18 

0. 1213 

7 

0.0749 

19 

0.1222 

8 

0.0752 

20 

0. 1248 

9 

0.0800 

21 

0. 1280 

10 

0 . 0862 

22 

0.1337 

1  1 

0.0879 

23 

0. 1419 

12 

0.0970 

24 

0.1440 

25 

0.1452 

37 

0 . 2567 

26 

0.1472 

38 

0 . 2578 

27 

0.1548 

39 

0 . 2597 

28 

0.1607 

40 

0.2730 

29 

0.1607 

41 

0 . 2798 

30 

0.1633 

42 

0.3443 

31 

0.1695 

43 

0 . 3493 

32 

0.1707 

44 

0 . 3493 

33 

0.1747 

45 

0.3607 

34 

0.1897 

46 

0 . 3967 

35 

0.2133 

36 

0.2162 

38 


TABLE  19. 

Ranks  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Strands 


Stress:  410.7  ksi 
Temperature:  110  deg  C. 
Number  of  specimens:  78 


1 

0.0095 

21 

0.1707 

41 

0.2487 

61 

0 . 3577 

2 

0.0113 

22 

0.1780 

42 

0.2513 

62 

0 . 3633 

3 

0.0217 

23 

0. 1825 

43 

0.2630 

63 

0 . 3663 

4 

0.0602 

24 

0.1848 

44 

0.2635 

64 

0 . 3673 

5 

0.0797 

25 

0. 1865 

45 

0.2650 

65 

0 . 3753 

6 

0.0855 

26 

0. 1880 

46 

0.2658 

66 

0 . 3758 

7 

0 . 0985 

27 

0.1915 

47 

0.2735 

67 

0.3802 

8 

0.0995 

28 

0.1958 

48 

0.2762 

68 

0 . 3858 

9 

0.1063 

29 

0.2013 

49 

0.2828 

69 

0 . 3942 

10 

0.1232 

30 

0.2032 

50 

0.2928 

70 

0 . 3977 

1  1 

0.1272 

31 

0.2083 

51 

0.2947 

71 

0.4268 

12 

0.1315 

32 

0.2170 

52 

0.2957 

72 

0.4433 

13 

0.1347 

33 

0.2195 

53 

0.3008 

73 

0.4478 

14 

0.1423 

34 

0.2217 

54 

0.3068 

74 

0.4745 

15 

0.1432 

35 

0.2235 

55 

0.3087 

75 

0.5010 

16 

0. 1453 

36 

0.2313 

56 

0.3193 

76 

0.51 15 

17 

0. 1497 

37 

0.2345 

57 

0.3252 

77 

0.5518 

18 

0. 1603 

38 

0 . 2437 

58 

0.3267 

78 

0 . 5685 

19 

0. 1665 

39 

0 . 2463 

59 

0.3278 

20 

O 

CD 

00 

ru 

40 

0.2478 

60 

0.3492 

39 


TABLE  20. 

Ranks  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Strands 

Stress:  410.7  ksi 
Temperature:  100  deg  C. 

Number  of  specimens:  53 


1 

0.039a 

15 

0.2882 

29 

0.5972 

43 

1.1030 

2 

0.0815 

16 

0 . 2983 

30 

0.6042 

44 

1 . 1370 

3 

0.0875 

17 

0.2997 

31 

0.6465 

45 

1.1380 

4 

0.1233 

18 

0.3087 

32 

0.7013 

46 

1 . 1750 

5 

0.1342 

19 

0.3183 

33 

0.7042 

47 

1 .2720 

6 

0.1787 

20 

0 . 3697 

34 

0.7808 

48 

1.3200 

7 

0. 1833 

21 

0.4415 

35 

0.8052 

49 

1 .4800 

8 

0.1880 

22 

0.4620 

36 

0.8473 

50 

1 .8320 

9 

0.1982 

23 

0 . 4625 

37 

0.8567 

51 

2.0690 

10 

0.2533 

24 

0.5312 

38 

0.8732 

52 

2.0700 

1  1 

0.2697 

25 

0.5535 

39 

0.8743 

53 

2. 1640 

12 

0.2697 

26 

0.5550 

40 

0.9833 

13 

0.2782 

27 

0.5550 

41 

0.9907 

14 

0.2783 

28 

0 . 5837 

42 

1 .0070 

TABLE  HI . 

Ranks  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Slrands 

Stress:  373.9  ksl 
Temperature:  120  deg  C. 

Number  of  specimens:  80 


1 

0. 1427 

21 

0.6520 

41 

1 .2043 

61 

1 .5947 

2 

0.1782 

22 

0.6600 

42 

1 .2063 

62 

1 . 6583 

3 

0.2633 

23 

0 . 6762 

43 

1 .2092 

63 

1 . 6828 

4 

0 . 2967 

24 

0.7008 

44 

1 . 2247 

64 

1 . 6887 

5 

0.3900 

25 

0.7923 

45 

1 .2278 

65 

1 .8027 

6 

0.H217 

26 

0.81 12 

46 

1.2411 

66 

1 .8480 

7 

0.4278 

27 

■  0.8143 

47 

1 .2800 

67 

1.9000 

8 

0.4318 

28 

0 . 8268 

48 

1 .3070 

68 

2. 1300 

9 

0.4375 

29 

0 .8585 

49 

1 .3091 

69 

2. 1980 

10 

0.4390 

30 

0.9097 

50 

1 .3345 

70 

2.2075 

1  1 

0.4542 

31 

0.9142 

51 

1 .3732 

71 

o 

OJ 

OJ 

12 

0.5313 

32 

0 . 9362 

52 

1.3913 

72 

2.3642 

13 

0.5330 

33 

0 . 9562 

53 

1 .4313 

73 

2.4338 

14 

0.5353 

34 

0 . 9566 

54 

1 .4352 

74 

2 . 6652 

15 

0.5397 

35 

0 . 9876 

55 

1.4415 

75 

2.8480 

16 

0.5538 

36 

0 . 9920 

56 

1 .4430 

76 

2 . 9662 

17 

0.5748 

37 

0.9990 

57 

1 .4980 

77 

3. 1533 

18 

0.5942 

38 

1 .0075 

58 

1  .5045 

78 

3.3112 

19 

0.6197 

39 

1.1120 

59 

1 .5045 

79 

3.7097 

20 

0.6323 

40 

1 . 1547 

60 

1 .5180 

80 

4.2721 

41 


TABLE  22. 

Ranl^s  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Strands 

Stress;  373.9  ksi 
Temperature:  110  deg  C. 

Number  of  specimens:  76 


1  0.0251 

2  0 . 0886 

3  0.0891 

4  0.2501 

5  0.3113 

6  0.3451 

7  0.4763 

8  0.5650 

9  0.5671 

10  0 . 6566 

11  0.6748 

12  0.6751 

13  0.6753 

14  0.7696 

15  0.8375 

16  0.8391 

17  0.8425 

18  0.8645 

19  0.8851 


20  0.9113 

21  0.9120 

22  0.9836 

23  1.0483 

24  1.0596 

25  1.0773 

86  1.1733 

27  1.2570 

28  1 . 2766 

29  1.2985 

30  1.3211 

31  1.3503 

32  1.3551 

33  1.4595 

34  1.4880 

35  1 . 5728 

36  1.5733 

37  1.7083 

38  1.7263 


39  1.7460 

40  1.7630 

41  1.7746 

42  1.8275 

43  1.8375 

44  1.8503 

45  1.8808 

46  1.8878 

47  1.8881 

48  1.9316 

49  1 . 9558 

50  2.0048 

51  2.0408 

52  2.0903 

53  2.1093 

54  2.1330 

55  2.2100 

56  2.2460 

57  2 . 2878 


58 

2.3203 

59 

2.3470 

60 

2.3513 

61 

2.4951 

62 

2.5260 

63 

2.9911 

64 

3.0256 

65 

3.2678 

66 

3.4045 

67 

3.4846 

68 

3.7433 

69 

3.7455 

70 

3.9143 

71 

4.8073 

72 

5.4005 

73 

5.4435 

74 

5.5295 

75 

6.5541 

76 

9.0960 

42 


TABLE  B3. 


Ranks  and  Breaking 

Times  ( 

Stress:  373, 

.9  k  s  i 

Temperature ; 

;  100  de 

g  c. 

Numb 

er  of  specimens 

:  76 

1 

0.  1050 

20 

1  .7042 

2 

0.1367 

21 

1.7267 

3 

0. 1516 

22 

1 .7276 

4 

0. 2833 

23 

1  .8098 

5 

0.3013 

24 

1 .8292 

6 

0 . 3882 

25 

2.0078 

7 

0.3902 

26 

2.2367 

8 

0.4083 

27 

2.7045 

9 

0.4378 

28 

2.7721 

10 

0.5113 

29 

2.8795 

1  1 

0.5863 

30 

2.9480 

12 

1  .  1243 

31 

2.9503 

13 

1.2750 

32 

3.0147 

14 

1  .2900 

33 

3.0525 

15 

1  .4228 

34 

3.  1625 

16 

1  .4651 

35 

3.3065 

17 

1.4882 

36 

3.3260 

18 

1 .5178 

37 

3.4085 

19 

1 .6163 

38 

3.4562 

)  for  Kevlar  H9/Epoxy  Strands 


39 

3 . 6768 

58 

6 . 6752 

40 

3.7078 

59 

6.7150 

41 

3.7095 

60 

6.8033 

42 

3.7997 

61 

6.8560 

43 

4.0242 

62 

6.8616 

44 

4.0253 

63 

7.2442 

45 

4.2678 

64 

7.5138 

46 

4.3113 

65 

8.0610 

47 

4.4470 

66 

8.3400 

48 

4.4788 

67 

8.4076 

49 

4.5420 

68 

8.4940 

50 

4.6260 

69 

9.0733 

51 

4.9577 

70 

9.4785 

52 

5.2326 

71 

10.2180 

53 

5.9757 

72 

1  1  .2280 

54 

6.0613 

73 

1 4 . 7858 

5^ 

6.0955 

74 

17.2493 

56 

6.1063 

75 

19.8447 

57 

6.6546 

76 

23.0297 

43 


TABLE  24. 

Ranks  and  Breaking  Times  (hr)  for  Kevlar  49/Epoxy  Strands 


Stress:  343.2  ksi 
Temperature:  120  deg  C. 
Number  of  specimens:  73 


1 

0.7367 

20 

3.9653 

2 

1.1627 

21 

4.2488 

3 

1 . 8945 

22 

4.3017 

4 

1 .9340 

23 

4 . 3942 

5 

2.3180 

24 

4.6416 

6 

2.6483 

25 

4.7070 

7 

2.8573 

26 

4.8885 

8 

2.9918 

27 

5. 1746 

9 

3.0797 

28 

5 . 4962 

10 

3.1152 

29 

5.5310 

1  1 

3. 1335 

30 

5 . 5588 

12 

3.2647 

31 

5.6333 

13 

3.4873 

32 

5.7006 

m 

3.5380 

33 

5.8730 

15 

3.6335 

34 

5.8737 

16 

3 . 654 1 

35 

5.9378 

17 

3.7645 

36 

6.1960 

18 

3.8196 

37 

6.2217 

19 

3.8520 

38 

6.2630 

39 

6.3163 

58 

9.9316 

40 

6.4513 

59 

10.0180 

41 

6.8320 

60 

10.4028 

42 

6.9447 

61 

10.4188 

43 

7.2595 

62 

10.7250 

44 

7.3183 

63 

10.9411 

45 

7.3313 

64 

1  1 .7962 

46 

7.7587 

65 

12.0750 

47 

8.0393 

66 

12.6933 

48 

8.0693 

67 

13.5307 

49 

8. 1928 

68 

13.8105 

50 

8.4166  ■ 

69 

14.5067 

51 

8.7558 

70 

15.3013 

52 

8.8398 

71 

16.2742 

53 

9.2497 

72 

18.2682 

54 

9.2563 

73 

19.2033 

55 

9.5418 

56 

9.6472 

57 

9.6902 

44 


TABLE  25. 

Ranks  and  Breaking  Times  (hr) 

Stress:  3H3 .2  k  s I 
Temperature:  110  deg  C. 

Number  of  specimens:  54 


1 

3.4347 

15 

10.9848 

2 

5.0883 

16 

1  1  .  1485 

3 

5.2947 

17 

1 1 .4787 

4 

5.8670 

18 

1 1 .4930 

5 

7.1213 

•  19 

1 1 .6778 

6 

8.5342 

20 

1 2 . 6927 

7 

8.6980 

21 

12.9188 

8 

8.7880 

22 

13.4178 

9 

9.2728 

23 

14.3510 

10 

9.6667 

24 

14.8937 

1  1 

9.8350 

25 

15.0912 

12 

9.9262 

26 

15. 1522 

13 

10.0907 

27 

16.6075 

14 

10.6818 

28 

17.5983 

for  Kevlar  49/Epoxy  Strands 


29 

18.0317 

43 

30.5233 

30 

19.2267 

44 

30.5433 

31 

19.4867 

45 

30.6183 

32 

19.5367 

46 

30.9917 

33 

19.7300 

47 

31 .4000 

34 

21 .5833 

48 

35.0900 

35 

22.5367 

49 

37.8417 

36 

22.8000 

50 

39.2433 

37 

24.4950 

51 

40.8633 

38 

26.4000 

52 

45.  1067 

39 

27.0783 

53 

60. 1983 

40 

28.1317 

54 

63.6067 

41 

28.5900 

42 

28.6517 

45 


TABLE  B6. 

Ranks  and  Breaking  Times  (hr)  for  Kevlar 

Stress:  343. B  ksi 
Temperature :  100  deg  C. 

Number  of  specimens:  58 


1  4.5145 

2  7.4662 

3  11.3375 

4  12.9017 

5  14.5092 

6  18.2700 

7  18.2867 

8  18.4283 

9  18.9767 

10  19.6067 

11  21.0517 

12  24.5867 

13  25.2250 

14  25.4050 

15  26.1517 


16  26.1750 

17  26 . 9767 

18  27.2833 

19  27.8417 

20  27.9033 

21  28 . 5850 

22  29 . 7867 

23  33.2800 

24  35.0933 

25  35.3183 

26  36.7000 

27  36.9050 

28  39.6950 

29  41.1783 

30  41.4117 


31  41.6983 

32  41.9100 

33  42.0117 

34  43.9867 

35  46.0433 

36  47.1717 

37  47.6983 

38  48.6550 

39  49.3833 

40  51.9239 

41  52.2283 

42  54.1567 

43  57.8333 

44  59.4983 

45  60.6517 


49/Epoxy  Strands 


46  65.8350 

47  65.9333 

48  70.9167 

49  74.5300 

50  74.7133 

51  87.0617 

52  102.2000 

53  103.5917 

54  105.4233 

55  109.3800 

56  138.4283 

57  156.8033 

58  168.6850 


46 


TABLE  E7. 

Ranks  and  Breaking  Times  (hr) 

Stress:  BBS. I  ksi 
Temperature:  1 10  deg  C. 

Number  of  specimens:  59 


1 

58. 

16 

164 

E 

66. 

17 

168 

3 

CD 

CD 

18 

168 

4 

IE. 

19 

168 

5 

74. 

BO 

17B 

6 

CD 

cn 

B1 

186 

7 

lOB. 

EE 

19B 

8 

110. 

B3 

196 

9 

o 

ru 

B4 

198 

10 

IBO. 

B5 

B09 

1 1 

IBO. 

B6 

B16. 

IB 

13B. 

El 

B3B, 

13 

135. 

B8 

B36. 

14 

143. 

B9 

B37. 

15 

144. 

30 

B40, 

for  Kevlar  49/Epoxy  Strands 


31 

B40. 

46 

360 

3B 

B5B. 

47 

360 

33 

B57. 

48 

360 

34 

B66. 

49 

360 

35 

B75. 

50 

376 

36 

B78. 

51 

384 

37 

B85. 

5B 

40B 

38 

B9B. 

53 

408 

39 

3B7. 

54 

418 

40 

336. 

55 

43B 

41 

336. 

56 

46B 

4B 

337. 

57 

463 

43 

338. 

58 

48B, 

44 

338. 

59 

5B3 

45 

357. 

47 


The  corresponding  estimates  of  the  vaV'iances  of  the  estimators  a(a,T)  and 

^  /V  ^ 

b(a,T),  and  of  the  covariance  between  a{a,T)  and  b(a,T),  are 
^  A  X  3  3  ^  ^  9  9 

Var(a(a,T))=  I  I  c^.c^.s^^,  Var(b(a,T))  =  I  I  c^c^s-j,  and 

^  ^  a''=’>  3=1  3  9  i=4  j=4 

Cov(a(a,T),  b(a,T))  =  ^i^j^ij»  ^2*  ‘^4»  ^5’  <^6»  ^7*  ^8* 

-1  2  -1  -2  -1 

Cg)  =  (  1,  a,  T  ,  1,  a,  a  ,  T  >  t  ,  ot  ).  Subsequent  formulas  for 

^  A 

estimates  of  the  variances  of  yJa,T)  and  R^.(a,T),  and  for  associated 

confidence  sets  for  y  (o.t)  and  R.(a,T),  are  identical  to  those  presented  in 

P  ^ 

Appendix  A,  with  (a,T)  substituted  for  a  throughout. 


G.  Appendix  C 

Results  analogous  to  those  found  in  Appendices  A  and  B  are  presented  here 
to  facilitate  computation  of  quantities  introduced  in  part  C. 

Table  28  gives  the  burst  pressure  data  for  the  seven  spools  used  in  this 
study.  (Spool  7  here  refers  to  the  old  spool  8.  The  original  spool  7  data 
are  omitted.) 

Tables  29  through  33  give  the  vessel  lifetime  data. 


2o 


The  natural  logarithm  of  the  ML  estimate  tp(a)  is  expressible  as 

yp(cr)  =  ®9  ®3^1  ^  ®5^3  ^  ®6^4  *  ®7^5  ^  ®8^6^  ^  ®10*^  ^  ®2'*^p® 

The  joint  probability  distribution  of  the  ML  estimates  (e-j,  ...,  9^q)  is 

approximately  10-dimensional  multivariate  normal  with  mean  vector  (e-i*  .... 

0ig)  and  covariance  matrix  z  estimated  by  z  =  (s^j),  i,j  =  1>  ...»  10, 
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Table  28.  Burst  pressure,  in  MPa,  of  NASA  vessels. 


Spool  1 

Spool  2 

Spool  3 

Spool  4 

Spool  5 

Spool  6 

Spool 

33.09 

32.75 

32.47 

36.03 

31.99 

34.47 

34.13 

34.89 

33.72 

34.54 

37.23 

33.72 

34.75 

35.32 

33.85 

34.99 

34.96 

35.58 

33.96 

35.16 

35.68 

35.58 

36.82 

36.54 


TABLE  29. 


Failure  Times  (hr)  for  NASA  Vessels 
Pressure  :  4.28  hs i 


Temperature 

:  25  deg 

C. 

Number  of  specimens: 

:  31 

Spool  1 

Spool  2 

Spool  3 

Sp  0  0 

444.4 

2.2 

12.5 

254 

755.2 

8.5 

14.6 

1  148 

952.2 

9.  1 

18.7 

1569 

1108.2 

10.2 

101.0 

1750 

22.  1 

1802 

55.4 

111.4 

158.7 


4 

Spool  5 

Spool  6 

Spool  7 

1 

8.3 

CD 

CD 

00 

ru 

5 

13.3 

15.0 

590. 1 

3 

87.5 

144.0 

638.2 

6 

243.9 

1 
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TABLE  30. 

Failure  Times  (hr)  for  NASA  Vessels 
Pressure:  3.97  ksi 
Number  of  specimens:  24 


Spool  1 

Spool  2 

Spool  3 

453.4 

71  .2 

19.  1 

664.5 

199.1 

24.3 

930.4 

403.7 

69.8 

1755.5 

453.4 

136.0 

514. 1 

544.9 

694.1 

Spool  4  Spool  5  Spool  6 

876.7  541.6  514.2 

1275.6  1254.9 

1536.8 
6177.5 


Spool  7 


554.2 

2046.2 
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TABLE  31 . 

Failure  Times  (hr)  for  NASA  Vessels 
Pressure:  3.68  ksi 
Temperature:  25  deg  C. 

Number  of  specimens:  23 


Spool  1 

Spool  2 

Spool  3 

Spool  4  Spool  5 

Spool  6 

1 1487.3 

1134.3 

1087.7 

13501.3  11727.1 

225.2 

14032. 

1824.3 

2442.5 

29808. 

6271 . 1 

30977. 

1920.1 

7996.0 

2383 . 0 

3708.9 

5556 . 0 

g  Grouped 

time  in  the 

i  n  t  e  r  V  a  1 

(4872. ,4944. ) . 

Spool  7 

2974.6 
9 

7332 . 0 

7918.7 

9240.3 

9973.3 


TABLE  32. 

Failure  Times  (hr)  for  NASA  Vessels 
Pressure:  3.38  l^si 
Temperature:  25  deg  0;,. 


Number  of 

specimens 

:  El 

Spool  1 

Spool  2 

Spool  3 

spool  H 

Spool  5 

Spool  6 

Spool  7 

14440. 

8615. 

9120. 

7320. 

20231 . 

16104. 

35880. 

20233. 

Censored 

Times  and  Quantities 

Spool  1 

Spool  2 

Spool  3 

Spool  4 

Spool  5 

Spool  6 

Spool  7 

H4376(4) 

3264 ( 1 ) 

19E(  1 ) 

44376(3) 

44376(4) 
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TABLE  33. 

Censored  Times  (hr)  and  Quantities  for  NASA  Vessels 

PressuretB.'+Bksi 

Temperature :  25  deg  C . 

Number  of  specimens:  40 


Spool  1 

44544(2) 

44688(3) 


Spool  2 

44544 ( 1 ) 
44688(5) 


Spool  3 
44688(3) 


Spool  4 
44688(7) 


Spool  5 

44544  (  1  ) 
44688(7) 


Spool  6 
44688(7) 


Spool  7 

1416(1) 

44688(3) 
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H.  Figures 

Figures  1  through  19  display  many  of  the  statistical  inferences  presented 
in  parts  A,  B,  C,  and  D. 


-56- 


Fiber  Stress 


Maximum  likelihood  est 
failure  probabilities 
90  (Room  temperature,  I 


data  on  Kevlar/ epoxy 
tes  for  percentile  fa 
0,  and  50  (room  tempe 


Max 

for 


3102 


estimates  of  lifetimes  for  Kevlar/epoxy  strands  under  tensi 


Figure  9.  ML  estimates  of  lifetimes  for  Kevlar/epoxy  strands  under  tensi 


Figure  13.  ML  estimates  of  lifetimes  for  Kevlar/epoxy  strands  under  tens 


verage 


ML  estimates  of  lifetimes  f 
based  on  data  from  high  stn 
law  fit.  (Room  temperature 


O-J-,  and  .50  failure  probabilities,  bas 
and  power  law  fit. 
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